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Abstract. Transport properties of the charge ordering compound /3— Nao.33V20s are studied in the tem- 
perature range from 30 K to 300 K using current driven DC conductivity experiments. It is found that 
below the metal-insulator transition temperature (Tmi = 136 K) this material shows a nonlinear charge 
density modulation behavior. The observed conductivity is discussed in terms of a classical domain model 
for charge density modulation transport. 

PACS. 71.30.+h Metal-insulator transitions and other electronic transitions - 72.20.Ht High-field and 
nonlinear effects - 72.80.-r Conductivity of specific materials 



e-mail: 



1 Introduction 

There are many solid state materials which are one - di- 
mensional and metallic at room temperature. Due to the 
coupling of the electrons to the underlying lattice the 
metallic state in these materials is usually not stable, lead- 
ing to a phase transition into a charge modulated state at 

. low temperatures. The low temperature ground state of 
the coupled electron-phonon system is characterized by 

, a gap in the single-particle excitation spectrum, by col- 
lective modes formed by the electron-hole pairs, and by 
a deformation of the lattice. Many of those compounds 
are inorganic (K . 3 MoO 3 , NbSe 3 , TaS 3 ) 0EE], but these 
properties are found in organic compounds as well 
((2.5(OCH 3 ) 2 DCNQI) 2 Li, TTF-TCNQ) gHHj. The best 
known type of this class of phase transitions is the Peierls 
transition |B] . In this case the material develops the charge 
density modulation state at low temperature with a si- 
multaneous deformation of the lattice, and the opening of 

, a gap in the quasiparticle excitation spectrum. As long 
as the charge density modulation is pinned by lattice or 
impurities, these materials can be described as narrow- 
band-gap semiconductors or even insulators. An excep- 
tion to this rule is NbSe 3 , which remains a semimetal at 
low temperature [7] . One of the more intriguing features of 
low dimensional charge ordered materials are their nonlin- 
ear transport properties. In the charge density modulation 
ground state they only exhibit ohmic conductivity below 
a certain threshold applied electric field Et- Below this 
field, the charge density modulation is pinned by lattice, 
impurities, defects, and grain boundaries and the conduc- 
tivity is solely due to a strongly temperature dependent 
quasiparticle transport. For fields above Et, the conduc- 
tivity becomes strongly enhanced and nonlinear due to 



the contribution of the now moving depinned charge den- 
sity modulation 8.9 l7)|lll) . In addition to the nonlinear 
behavior, charge density modulated materials often show 
an alternating current response to a static applied field. 
This latter may be either due to so called ratcheting of 
the charge density modulation phase as observed in for 
instance NbSe 3 jjj], generally referred to as narrow band 
noise, or due to macroscopic polarization oscillations as 
observed in for instance blue bronze at low temperatures 

The recently revived interest in the vanadium bronze 
j3— Nao. 33 V205 has been triggered by the observation of 
a charge ordering transition and sparked once more 
by the observation of the pressure induced superconduc- 
tivity |14j in this electronically low dimensional material. 
The vanadium bronze f3— Nao. 33 V20s has been the sub- 
ject of various structural studies during the last 40 years 
|15II16| . At room temperature /?— Nao. 33 V205 is a highly 
anisotropic metal with a site occupancy disorder of the 
sodium atoms. Around Tjy a = 240 K a second order phase 
transition occurs leading to an ordering of the sodium 
atoms and a doubling of the primitive cell along the b di- 
rection |17| . A charge ordering transition, which is thought 
to be driven by the electron phonon coupling, occurs at 
Tmi = 136 K and is accompanied by a further tripling 
of the unit cell along the b direction |17U18j . leading to a 
commensurate charge modulated state with a period of 6b. 
Temperature dependent measurements of the magnetic 
susceptibility revealed a magnetic transition from a para- 
magnetic to an antifcrromagnetic state at Taf — 22 K 
|13| . Finally, as mentioned above, (3— Nao. 33 V20s becomes 
a superconductor for pressures above 8 GPa and low tem- 
peratures (8 K) as revealed by recent pressure dependent 
resistivity measurements |14) . Optical conductivity data 
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suggest that this system may be understood in terms of a 
small polaron model, where the charge ordering in fact cor- 
responds to an ordering of the polarons |I9II20| . Since pre- 
vious experiments strongly suggest that {3— Nao.33V20s 
is a low dimensional conductor with an electron-phonon 
interaction induced charge ordering, this material should 
also show the usual nonlinear transport properties dis- 
cussed above. Indeed the results presented in this paper of 
the first detailed measurements of the field dependence of 
conduction in the sodium bronze (3— Nao.33V20s are fully 
consistent with this picture. The observed nonlinear trans- 
port properties are well described using a classical domain 
model. The observed charge density modulation conduc- 
tivity increases with increasing temperature, suggesting a 
screening of the charge density modulation pinning by the 
thermally excited carriers. 



2 Temperature dependent transport 

Single crystal samples have been prepared as described 
elsewhere jT^J . Platelets with typical dimensions 4x2 
x 0.2 mm 3 were mounted on the cold finger of a flow 
cryostat and contacted using 50 /im diameter platinum 
wires. Four wires were fixed on the sample surface using 
silver paste, spaced 1 mm apart. Measurements were per- 
formed along the b axis in a four probe configuration us- 
ing a Keithley 236 source-meter in a current driven mode. 
The resistivity was calculated from the measured resis- 
tance and the geometry of the sample. 

First we turn to the temperature dependence of the 
low-field ohmic resistance, of which a typical example is 
shown in FigQ] 31] • The strong increase of the resistiv- 
ity below Tmi — 136 K clearly signals the charge order- 
ing transition, consistent with results published in liter- 
ature. Note that also the sodium ordering transition at 
Tno, = 240 K leads to a change in the temperature de- 
pendence of the resistivity (see also inset Fig. QJ). The 
small enhancement of the conductivity below T^ a can be 
understood in terms of the decreased amplitude of the spa- 
tial potential variations on the vanadium sites due to the 
ordering of the sodium subsystem |TJj]. Although the resis- 
tivity at high temperatures (T > Tmi) is fairly small, as 
expected for a bad metal, it does not exhibit the metal- 
lic behavior as observed previously |17|. The reason for 
this could be that the samples are slightly misaligned. In 
this case the measured resistivity is a combination of the 
resistivity of the metallic b-axis and the insulating per- 
pendicular axes, leading to the observed temperature de- 
pendence. X-ray diffraction experiments on the samples, 
however, have ruled out this option in showing that to 
within a degree the samples are indeed b-oriented. A more 
likely reason is that the sodium stoichiometry of the sam- 
ple slightly deviates from x = 0.33, which is known to lead 
to a rapid loss of metallic behavior and eventually to the 
disappearance of the metal-insulator transition |I3| . Such 
deviations lead to an additional disorder in the sodium 
site occupancy, which in turn leads to a more disordered 
potential on the vanadium sites which make up the one- 
dimensional conduction chains |19] . Since in particular low 
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Fig. 1. Resistivity as a function of temperature in 
p— Nao.33V205. Inset: change in conductivity near the sodium 
ordering transition at Tjv a ~ 240 K. 



dimensional systems are very susceptible to disorder, this 
may lead to a small disorder induced gap and to the ob- 
served non-metallic behavior. 

An estimate for the transport gaps above and below 
the charge ordering transition temperature is obtained fit- 
ting a simple activated behavior, p — poexp(AfkT), to 
the data. Well within the charge ordered phase (30 K < 
T < 80 K) we find po = 630 mJ7cm and a gap Alt — 
548 K. At high temperatures (150 K < T < 300 K) we find 
po = 75 mJ?cm and a gap Ajjt — 472 K. Surprisingly, the 
transport gap in the charge ordered phase is only about 
15% larger than the one found for the high temperature 
phase. The major change is found in the prefactors which 
differ by an order of magnitude. This is consistent with 
the expectation that the majority of the charge carriers 
is frozen out by the charge ordering at Tmi- In a similar 
analysis Yamada et al. ^3] found an activation energy of 
538 K, in good agreement with the present value. In previ- 
ous work we found an optical gap of 2Z\=2450 K [TJ]. Most 
likely, this optical gap corresponds to the energy needed to 
free a quasiparticle from the charge ordered state. Clearly 
then, the transport gap must have a different origin. The 
origin for this difference could be that the observed trans- 
port gap is in fact a sodium disorder induced pseudo gap 
in the charge carrier spectrum, where the charge carriers 
themselves are present due to an incomplete charge order- 
ing, rather then due to thermal excitation as is usually the 
case. This would be consistent with the small difference 
in the transport gap values below and above the charge 
ordering transition, as well as with the observed finite low 
frequency conductivity in the optical data |19| . 

However, closer inspection of the low temperature con- 
ductivity shows that it does not strictly follow an activated 
behavior. In particular, below 60 K there seems to be an 
enhancement of the conductivity, which presumably is due 
to the presence of an additional conductivity channel. In a 
purely one-dimensional model, such an enhancement may 
originate from the presence of mid-gap bound states of 
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Fig. 2. Temperature dependent conductivity (symbols, same 
data as in Fig. 1.) together with a fit to an activated two par- 
ticle behavior, Eq. (1) (dark line). The gray line shows a fit 
of thermally activated quasiparticle transport to the high tem- 
perature part of the data. 



amplitude 7r-solitons and quasiparticles, as described by 
Brazovskii [221 • Indeed, the low temperature conductivity 
is better described by the empirical from 



a(T) = a° QP cxp(-A QP /kT) + a° s exp(-A s /kT). (1) 

Here the first term on the right hand side accounts for 
the contribution of the quasiparticles, whereas the second 
term takes the midgap state conductivity into account. 
In the one-dimensional description, the midgap states are 
located halfway the quasiparticle gap, i.e. As ~ Aqp/2. 
Fitting Eq. JTJ to the data yields A QP = 995 K; A s = 
487 K, and a ratio cr^ p /a° s ~ 0.4 x 10 3 . Note that the 
quasiparticle gap obtained in this way becomes compa- 
rable to the gap found in optical experiments [JJ]. The 
conductivity ratio shows that just below the phase tran- 
sition quasiparticle transport dominates the conductiv- 
ity, whereas the 'midgap' contribution becomes impor- 
tant at lower temperatures only, despite its smaller en- 
ergy gap. Finally, we note that a similar analysis to blue 
bronze, K0.3M0O3, data leads to similar conclusions. For 
both these cases, one can worry whether a one-dimensional 
model is applicable at all. Indeed, the midgap states de- 
scribed by Brazovski require the existence of topological w- 
solitons which do not exist in three dimensions. Therefore, 
the above analysis merely shows the presence of additional 
excitations halfway the gap. Though in sodium vanadate, 
these might originate from the sodium disorder, one would 
not expect a similar disorder in blue bronze. Therefore, the 
precise nature of the midgap states remains unresolved at 
present. One interesting thought is that they might result 
from excitations inside domain walls which separate the 
charge density modulation ordered regions in the samples. 



3 Field dependent transport 

The charge density modulation in low dimensional sys- 
tems is usually pinned by the underlying lattice, by im- 
purities or by structural defects. When the charge den- 
sity modulation is incommensurate with the underlying 
lattice, the main pinning centers are the impurities and 
other defects. In this case, the pinning might be consid- 
ered relatively weak, and the charge density modulation 
can fairly easily move in the material upon application of 
an external electric field, leading to the typical non-linear 
charge density wave conduction. When the modulation is 
commensurate with the lattice, the pinning is considered 
to be strong, usually much stronger than the impurity pin- 
ning, and the large pinning energy prevents conduction of 
the charge density modulation. The charge modulation in 
j3— N&o. 33V 2O 5 is commensurate with the underlying lat- 
tice. Since the modulation period is, however, rather large 
(6b) the commensurability pinning in this system is ex- 
pected to be relatively weak, opening the possibility of 
non-linear charge density wave conduction in the charge 
ordered phase. 

To study the nonlinear transport properties of single 
crystal j3— Nao.33V2 0s, current driven field dependence 
measurements |23| were performed in the temperature range 
65 K - 300 K. The obtained transport properties are very 
similar with those obtained in CDW systems. A typical 
example of the conductivity measured along the b-axis at 
65 K is displayed in Fig. 
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Fig. 3. Nonlinear field dependent conductivity in 
13— Nao.33V205 measured along the b axis at 65 K dis- 
playing three charge density modulation transport regimes. 

The behavior of the conductivity shows three regimes. 
Below 0.06 mV/cm (first threshold field) the conductiv- 
ity is field-independent and mainly due to quasiparticle 
transport. At about 0.06 mV/cm the conductivity shows 
a nonlinear increase resulting from an incoherent contribu- 
tion of the charge density wave to the conductivity. This 
behavior saturates around 4 mV/cm. Above the second 
threshold field at 30 mV/cm a steep increase of the con- 
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ductivity takes place signaling the onset of the coherent 
charge density modulation regime. Despite the strong in- 
crease of the conductivity, it never reaches values compa- 
rable to the metallic state. 

The values for the threshold fields are smaller than 
those found, for example in blue bronze. Zawilski et al. 
reported 40 mV/cm for the first threshold field and around 
2000 mV/cm for the second threshold field in K0.3M0O3 
at 60 K while Mihaly et al. [23| found 40 mV/cm in 
Ko^-^Na^MoOs around the same temperature. Fleming 
et al. |251 found the first threshold field around 500 mV/cm 
in TaS3 and 90 mV/cm in K0.3M0O3 both measured at 
60 K. In Ko.s-^Na^MoOg (x = 0, 0.02, 0.05, 0.1), at 77 K, 
Wang et al. EE] reported values of (300, 670, 750, 950) 
mV/cm for the first threshold field at 180 K. Beauchene et. 
all reported similar values in Rbo.sMoOV Kiintscher 
et al. [2Hj found in K0.3M0O3 and Rbo.3Mo03 values of 
150 mV/cm for the first threshold field at 60 K while van 
Loosdrecht et al. 10 found 200 mV/cm at the same tem- 
perature. Although the values of the first threshold field 
for K0.3M0O3 show some variation, they are typically in 
the 50-200 mV/cm range, i.e. substantially higher than 
the present values observed for j3— Nao.33V20V This is 
consistent with the notion that the temperature depen- 
dence of the low field transport observed in (3— Nao.33V205 
is due to incomplete charge ordering resulting from disor- 
der in the sodium sublattice. The presence of charge car- 
riers, even at low temperatures, leads to screening of the 
pinning potential, and hence to a lowering of the threshold 
field. In contrast, the presence of charge carriers in blue 
bronze is almost entirely due to thermal quasiparticle ex- 
citations. 
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Fig. 4. Nonlinear conductivity along the b axis as a function 
of electric field in f3— Nao.33V20s measured from 65 K to 300 K 
(in 5 K steps between 65 K - 125 K, 140 K - 230 K, and 250 K 
- 300 K; in 2 K steps in the transition temperatures regions: 
128 K - 146 K, and 236 K - 246 K). 

Figure 0] displays the electric field dependence of the 
conductivity along the b-axis for a variety of tempera- 
tures between 65 K and 300 K. The nonlinear behavior 



is most pronounced at low temperatures, and is slowly 
decreasing upon increasing the temperature towards the 
transition temperature Tmi- The first threshold field Et, 
i.e. the field required to induce charge density modulation 
conductivity, is observed below Tmi only. As the tempera- 
ture is reduced, the threshold fields increase, evidencing a 
strengthening of the charge density modulation pinning. 
This is due to the reduction of the free carrier concen- 
tration at lower temperatures, leading to a less effective 
screening of pinning centers. The second transport regime, 
the incoherent moving regime, shows an increasing con- 
ductivity followed by saturation upon increasing field. Go- 
ing up in temperature, the field at which this saturation 
is reached increases, until around 90 K it merges with the 
second threshold field. 

Even above the charge ordering transition temperature, 
the nonlinear behavior has not entirely disappeared, al- 
though there are no clear threshold fields anymore. Non- 
linear conduction is still observed up to the sodium or- 
dering temperature T^ a — 240 K, consistent with the dis- 
order induced non-metallic behavior observed above Tmi- 
Finally, above 240 K a nearly field independent conduc- 
tivity is observed. 
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Fig. 5. Temperature dependence of the threshold fields for 
T < Tmi, obtained from the conductivity data. Open symbols 
(left scale): first threshold field; Filled symbols (right scale): 
second threshold field. 

The temperature dependencies of the first and second 
threshold fields are shown in Fig. The open symbols 
display the first threshold field, Et, obtained from the 
conductivity curves by taking the values where the con- 
ductivity starts to increase. The lower part displays the 
second threshold field, E T , obtained in a similar manner. 
Although there is a factor of 10 3 difference between the 
two threshold fields, the displayed behavior is qualitatively 
the same: the threshold field strongly decreases with in- 
creasing temperature. With some exceptions ( for exam- 
ple K0.3M0O3 or TaS3 ) this type of behavior for in 
particular the first threshold field is common for many of 
the known charge density wave systems. At low tempera- 
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ture, when the number of quasiparticles is reduced, large 
electric fields may build up around pinning centers and 
the value of E T is essentially determined by the amount 
of disorder in the system . Going up in temperature, the 
thermally excited quasiparticles tend to homogenize the 
electric field inside the sample. This effect is a natural 
source for the redistribution of the driving fields inside 
the sample. 

A variety of models were proposed to describe the 
nonlinear conductivity observed in charge density mod- 
ulation materials, mostly based on the original suggestion 
of Frohlich (30] that conductivity is dominated by slid- 
ing CDW transport. Between them, there are two main 
approaches in describing the CDW conductivity [Dl$ll31| . 
The first one, treats the CDW as a classical particle which 
is moving in a periodic potential, with the period deter- 
mined by the period of the CDW [521155] . This model gives 
a sharp threshold field, Et for the onset of nonlinear- 
ity and a saturation of the conductivity for high electric 
fields. The second model, proposed by Bardeen and re- 
ferred to as the tunneling model [2], assumes that the 
nonlinear transport occurs as a result of coherent tun- 
neling of the CDW over macroscopic distances. Besides 
the threshold field, Et, the model gives another charac- 
teristic field Eg, which can be interpreted as a tunnel- 
ing barrier. The onset of the conductivity reveals a sharp 
threshold field and at high electric fields the conductivity 
saturates. Both models rely on a T = treatment of the 
problem, though nonzero temperature models based on 
thermally assisted flux creep have been discussed as well 
35 3j. The models described above fail to account for 
the charge density modulation conductivity observed in 
/3— Nao.33V205. Fukuyama, Lee and Rice |3*rjll3*7| discussed 
the effect of the pinning centers on the charge-density- 
waves dynamics. They focused on phase fluctuations and 
show that, when the material is characterized by weak pin- 
ning, the system can be though to break up into domains. 
Here, we integrate this notion in a simple phenomenolog- 
ical model describing the observed nonlinear transport in 
(3— Nao.33V205, by treating the sample as a collection of 
interconnected domains. Each domain is characterized by 
its own conductivity and threshold field so that the sam- 
ple can be considered as a collection of parallel and series 
nonlinear conduction paths. At low applied fields (smaller 
then Et), excited quasiparticles dominate the conductiv- 
ity of a single domain, leading to the strongly temperature 
dependent ohmic behavior. All charge density modulation 
domains are pinned by lattice and defects and will not 
start moving until the applied field exceeds a certain crit- 
ical field. At higher applied fields, some domains start to 
become depinned, leading to an additional charge density 
modulation contribution to the conductivity. The charge 
density modulation still cannot move as a whole due to 
the domain structure resulting from grain boundaries and 
strong pinning centers. The sample is now in the incoher- 
ent moving regime. Finally for applied fields exceeding a 
second critical field, the charge density modulation may 
move as a whole (or at least a percolation path exists be- 
tween the contacts), and the conductivity becomes com- 



pletely dominated by the moving charge density modula- 
tion transport. Within this model, the second critical field 
is then a percolation threshold and the transport regime 
above this field is dubbed coherent moving regime. 

The model sketched above can be made more quantita- 
tive by specifying a model for the charge density modula- 
tion conductivity itself. Here we take one of the simplest 
approaches, and describe the charge density modulation 
within a domain as a charged particle moving due to an 
applied field E in a viscous medium pQ . The charge density 
modulation conductivity is then given by. 

a cdm (E) = <Tc E ~ ET 9(E - E T ) , (2) 

where the Heaviside function 9(E — Et) assures that the 
equation is also valid for E < Et- In the incoherent mov- 
ing regime the conductivity can then be modeled as a 
collection of interconnected charge density modulation do- 
mains, shunted by the free carrier conductivity. In general, 
this is still a complex system which solution would require 
detailed knowledge of domain properties and their connec- 
tivity. Here we will take a more qualitative and statisti- 
cal approach, and model the charge density modulation 
system as a statistically large number of interconnected 
domains. The conductivity of this network of domains to- 
gether with the free carrier conductivity will then yield 
the total conductivity. The network can be considered as 
a collection of parallel and series conduction paths. For 
the series connected domains of a single conduction path, 
the conductivity will have the same form as Eq. (j5J), but 
with an effective threshold field Et = Y^i an d effec- 
tive conductivity a c — (^c) 1 ) • The total charge 
density modulation conductivity is then given by 



Vcdm(E) = 2_,^ 



E — El, , 
~^0(E 



E T ), 



(3) 



where j runs over the parallel conduction pathways, each 
having their own effective threshold field E T and effective 
conductivity a° c . If there are a large number of conduc- 
tion paths, the above summation can be replaced by an 
integration, from zero to infinity, weighted by the statisti- 
cal distribution of threshold fields and conductivities. For 
simplicity, we assume a lorentzian distribution of width 
7, centered on E T for the threshold fields, and a constant 
effective conductivity a J c — a c for each path. The total 
conductivity, including the free carrier contribution, ctoj 
can then be equated as 



o-(e) = Co + < 



2e(arctan(e) + arctan(e )) 



(e + e ) (tt + 2 arctan(e )) 



(4) 



with e = E T /j, e = (E - E T )/j. In the limit j/E T 
3> 1, which is valid for the present experiment, the total 
conductivity becomes: 



cr(e) = <t + o", 



2earctan(e) - ln(l + e 2 ) 



(5) 
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where e = E/'y 

In the derivation of Eq. (JSj quite a few assumptions 
have been made. The most important one is that domains 
are formed in the charge density modulation material, 
within which the charge density modulation transport is 
coherent. This is mainly based on the observation of the 
slow onset of the moving conductivity in /?— Nao.33V205 
(see Fig. |3| and Fig.^]) as well as in for instance K0.3M0O3 
The assumption of parallel paths of series resistors 
is pretty robust; allowing for connectivity between these 
paths would only lead to additional parallel pathways for 
conduction. Of course, there is nothing known on the sta- 
tistical distribution of domain properties. Taking a sym- 
metric lorentzian (or gaussian) distribution simply makes 
the integration over domains tractable. Taking an asym- 
metric distribution would be more physical. Also typical 
domain sizes are not known. X-ray diffraction measure- 
ments |38] . however, have shown that the width of the 
superstructure peaks originating from charge ordering is 
comparable to the sharpness of the fundamental peaks, 
thereby setting a lower limit of the domain sizes to ^100 
nm. Finally, we did not allow for a distribution of effective 
conductivities for the conduction paths. Numerical simula- 
tions taking lorentzian distributions a J c have shown, how- 
ever, that the shape of the field dependent conductivity 
does not strongly depend on this. 

Now return to the data in Fig.|3] We have fitted Eq. JSJ 
to the field dependent conductivities, measured at differ- 
ent temperatures between 65 K and 136 K. Note that the 
only free parameters in the fits are the width of the distri- 
bution 7, and the charge density modulation conductivity 
<T C , since the normal conductivity, ao, can be obtained di- 
rectly from the low field data. At low temperatures, a good 
approximation of er c can be made using the field depen- 
dent conductivity data taking the conductivity difference 
between the ohmic regime and the saturation regime. The 
only remaining fit parameter in this situation is the width 
of the lorentzian distribution, 7. The fits generally show 
a good agreement with the data, and some typical results 
of the fitting are shown in Fig. |SJ 

The temperature dependence of the threshold field dis- 
tribution width, 7, obtained from the fits is displayed in 
Fig.0 

At low temperature, the width increases with increas- 
ing temperature. It shows a pronounced peak at 80 K, 
after which it starts to decrease becoming very small in 
the region of the MI transition. The decrease with increas- 
ing temperature observed above 80 K is what one might 
intuitively expect; the enhanced screening will decrease 
the typical domain threshold fields, thereby decreasing 7. 
Apart from the decrease of the average pinning potential, 
there will also be changes in the statistical distribution as 
the temperature is lowered. At low temperatures one ex- 
pects that there will be a larger number of domains with 
a relatively small pinning potential, reducing the distri- 
bution width. Therefore we believe that the increase of 
the width observed in the low temperature results from a 
competition of enhanced screening and the formation of 
larger, more strongly pinned, domains due to the coales- 
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Fig. 6. Fits of the domain model, Eq. Q (solid lines), to the 
nonlinear transport data (open symbols) at 65 K, 70 K and 
75 K. 
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Fig. 7. Temperature dependence of the width 7 of the 
lorentzian threshold field distribution obtained from fits of the 
data in Fig. 0] to the domain model, Eq. (jJJ . 



cence of small, weakly pinned, domains as the temperature 
increases. 

Figure |Hl shows the temperature dependence of the ef- 
fective charge density modulation conductivity, <j c , and 
the low field ohmic conductivity, ctq. The moving charge 
density modulation conductivity is found to be almost an 
order of magnitude larger than the quasiparticle contri- 
bution. The small kink in the charge density modulation 
contribution around 90 K, which probably again results 
from the competition between screening and coalescence. 

The temperature dependence of both contributions show 
a similar activated behavior. For the quasiparticles this 
has already been discussed in section [3 The usual inter- 
pretation of the activated behavior of the moving den- 
sity modulation contribution is in terms of thermally ac- 
tivated flux creep similar to the flux creep of the 
Abrikosov flux lattice in superconductors |39| . From fits 
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Fig. 8. Temperature dependence of the effective charge den- 
sity modulation conductivity, <f c (open symbols), and the low 
field ohmic conductivity <ro (filled symbols), obtained from fits 
of the data in Fig. 0] to the domain model, Eq. @ (see text). 



Fig. 9. Fits of the domain-percolation model, Eq. ||3J (solid 
lines) to the nonlinear transport data (open symbols) at 65 K, 
70 K and 75 K. The inset shows the behavior of the percolation 
threshold with temperature. 



to an activated behavior we estimate the low temperature 
(T < 90 K) activation energies for the quasiparticle and 
the moving density modulation contributions to be 740 K 
and 805 K, respectively. For the ohmic contribution, this 
is in good agreement with the earlier results (Sec.EJ). 

The above discussed model gives a phenomenologi- 
cal understanding of the ohmic and incoherent transport 
regimes. It does not, however, describe the coherent trans- 
port regime, where a sharp rise in conductivity takes place. 
The appearance of this second threshold field can be un- 
derstood as follows. As the field increases, more and more 
local pinning potentials will be overcome, increasing the 
number of domains which contribute to the conductivity. 
Eventually this will lead to the formation of a percolation 
path between the contacts. We thus propose that the sec- 
ond threshold is in fact a percolation threshold. For fields 
bigger than the second threshold field, E?p, the CDW will 
then coherently move along such percolation paths, lead- 
ing to the sharp rise in conductivity. This leads to an ad- 
ditional contribution to the conductivity of the form of 
Eq. (@J, so that the total conductivity now becomes 

„ 2earctan(e) — ln(l + e 2 ) 

o-(e) = cr + <?c — L + 

eir 

2e p (arctan(e p ) + arctan(e 0p )) — ^(tt-¥") 
(e P + eo P )(7T + 2 arctan(e 0c )) 

where e p = (E — E^)/j p , and 6q p — E^/j p . j p and 
<7 P are the width of the distribution and the percolation 
charge density modulation conductivity, respectively. We 
have fitted this last equation to the low temperature (T < 
90 K) field dependent conductivities and some typical re- 
sults of the fitting are shown in Fig. El 

Clearly the data follows Eq. © quite well. The tem- 
perature dependence of the upper threshold field (see inset 
Fig. El) closely follows the results presented in Fig. El The 



reason for the sharp decrease of the percolation threshold 
field upon increasing temperature is the same as for the 
incoherent moving threshold, as has been discussed in sec- 
tion El namely the increased screening of impurities upon 
increasing temperature. Finally, we note that the fitting 
is fairly insensitive to the distribution width as long as 
7 P << as expected for a percolation threshold. 



4 Conclusions 

We presented detailed nonlinear transport experiments on 
(3— Nao.33V205 in the temperature range 30 K - 300 K. 
The low field data in the charge ordered phase show that 
two types of excitations contribute to the transport. The 
charge density modulation quasiparticle gap is found to 
be about 700-800 K, and likely depends on the sodium 
stoichiometry. Evidence for a second type of excitation, 
with a gap of ~ 500 K, has been presented, although the 
exact origin of this excitation remains unclear at present. 
It might be either a bound state of collective charge den- 
sity modulation excitations like the phason, or possibly 
an excitation within the domain walls between ordered 
charge density modulation domains. A competing model 
for the thermally activated low field charge transport in a 
low dimensional system in the presence of disorder is the 
variable range hopping (VRH) model 40,41 . Analyzing 
the low temperature data in terms of VRH conductivity 
indeed leads to a reasonable agreement with this model as 
well, again with deviations at temperatures below 50 K. 
Since the present data can not distinguish between these 
models we adapted the most widely used model here as 
well. 

The field dependent data clearly show the charge den- 
sity modulation nature of the insulating phase, and is very 
similar to the transport in other well known charge density 
modulation materials like K0.3M0O3. The phenomenolog- 
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ical domain model for nonlinear transport in charge den- 
sity modulation materials presented here is found to be in 
good agreement with the experiments, and we believe that 
this model should be applicable to other semi-conducting 
charge ordered materials as well. The model could be im- 
proved by taking a more realistic model for the statistical 
distribution of domain properties, the single domain trans- 
port, and by allowing for field dependent domain proper- 
ties. In particular the pinning fields are expected to be field 
dependent, since at high moving velocities the charge den- 
sity modulation excitations may excite quasiparticle exci- 
tations resulting from charge density friction. Finally, we 
believe that in low free carrier density (i.e. semiconduct- 
ing) compounds like (3— Nao^A^Os and K0.3M0O3 much 
of the temperature dependence of the transport proper- 
ties, including the observed decrease of threshold fields 
upon raising temperature, can be understood in terms of 
enhanced screening of pinning centers by thermally ex- 
cited charge carriers. 
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